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Abstract 

It is shown that the translational degrees of freedom of a large variety of molecules, from light 
diatomic to heavy organic ones, can be cooled sympathetically and brought to rest (crystallized) 
in a linear Paul trap. The method relies on endowing the molecules with an appropriate positive 
charge, storage in a linear radiofrequency trap, and sympathetic cooling. Two well-known atomic 
coolant species, 9Be+ and ^''^'^Ba+, are sufficient for cooling the molecular mass range from 2 to 
20,000 amu. The large molecular charge required for simultaneous trapping of heavy molecules and 
of the coolant ions can easily be produced using electrospray ionization. Crystallized molecular 
ions offer vast opportunities for novel studies. 

PACS numbers: 42.50.Vk, 32.80.Pj, 33.80.Ps, 39.25. +k 
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After the enormous success achieved in the field of cold atom manipulation, significant 
efforts are under way to develop similar methods for molecules. Samples of trapped ultracold 
molecules could be used for performing studies of molecular structure, chemical reactions, 
quantum optics and molecular Bose-Einstein condensates. 

While methods for trapping of molecules, such as magnetic traps 0], electrostatic traps 
radiofrequency (Paul) or electromagnetic (Penning) traps or dipole traps jj] have been 
demonstrated and are in part already well developed, translational cooling of molecules is 
a field still under intense development. Direct laser cooling is not applicable due to lack of 
closed transitions. A technique demonstrated early on for coolingof neutral and charged 
molecules is by means of cryogenic buffer gas such as ^He, see e.g. Its extension to 

■^He Q allowed reaching temperatures below 1 K. Two recently developed methods are the 
)roduction of ultracold diatomic neutral molecules by photoassociation of ultra-cold atoms 
^, and the deceleration of polar molecules by time- dependent electric fields j^. 

Another powerful method is sympathetic cooling of " sample" particles of one species by 
an ensemble of directly cooled (often by laser cooling) particles of another species via their 
mutual interaction. This technique, first demonstrated for ions in Penning traps 

00, 

is being applied to an increasingly wide variety of neutral and charged particles (atoms, 
molecules, elementary particles) in various trap types for applications ranging from mass 
spectrometry to quantum computing. 

In both Penning and radiofrequency traps, first studies showed that molecular ions could 
be svmpatheticall y co oled (sc) by laser-cooled (Ic) atomic ions to temperatures in the range 
of several K Il3|, and mass ratios down to m^clm\r. = 2/3 were achieved lljl. An 



important aspect of ion traps is that for sufficiently strong cooling the formation of an ordered 
structure (Coulomb crystal) results 15, 3l- Recently it has been shown that sympathethic 
crystallization of molecular ions is possible in a linear rf trap where the molecules are stably 
incorporated into the atomic Coulomb crystal The molecular ions included MgH"*", 

MgD"*" (cooled by Mg'^) and ^^O^ (cooled by ^'^Ca^ or ^"^Mg^), i.e. mass ratios down to 
0.6. In sympathetic crystallization of atomic ions, a mass ratio range of 0.8 - 1.8 has been 
achieved |l9|,|2o|. 

The range of molecular mass that can be sympathetically cooled to temperatures of the 
order of 10 K, where the ion ensemble is still in a gas state, has very recently been debated 



in theoretical work 
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A molecular mass range from 8 to 192 amu (mass ratios 
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1/3 - 8) was found using molecular dynamics simulations to be accessible using ^^Mg"*" as 
coolant ion. 

In the present simulations we study a much larger mass range and the regime of much 
lower temperatures 1 K) and study whether sympathetic crystallization can be reached. 
We find that essentially all molecular masses can be sympathetically crystallized by one of 
two commonly used species of laser-coolable ions, ^Be"*" and ^^''Ba'*'. The only requirement 
for this general method is an appropriate charge state for the molecules. Single positive 
charges are sufficient for small molecules (mass 2 - 2000 amu). For heavier molecules higher 
positive charge states are required in order to allow for reliable simultaneous trapping. These 



are easily produced using electrospray ionization j2J]. The ability to store molecules in an 
almost motionless state in a collision-free ultra-high vacuum environment for essentially 
unlimited time is expected to open up vast opportunities for high precision spectroscopy 
and the study of slow molecular processes. 

In order to model sympathetic cooling in rf traps, it is crucial to take collisions into ac- 
count precisely. Noninteracting ions perform an oscillation at the frequency of the applied 
rf field, but their cycle-averaged energy is constant if ion-ion interactions can be neglected. 
When interactions, i.e. collisions, between the ions are taken into account, energy gain from 
the rf field (rf-heating) and transfer of energy from one species to another, i.e. sympathetic 



heating and coo 



mechanisms 
tia 







ing, occurs. Approximate models have been proposed to describe these 



22| . However, it is highly desirable to perform calculations that are essen- 
ly free of approximations. We therefore choose a molecular dynamics (MD) approach 

Mm- 

The simulations are based on solving Newton's equations of motion for the laser-cooled 
(Ic) and sympathetically cooled ions (sc) 

m^Xi = giEtrap(Xi, t) + Fc,i({Xj}) + FL(Xi, t) , (1) 

where i = 1, . . . , Ni^ + Nsc (A^k and Age are the numbers of laser-cooled and sympatheti- 
cally cooled particles, respectively). Positions, charges and masses are Xj, Qi, and mj, the 
Coulomb force Fc,j = {Qi/4:7ieo)ViJ2jQj/''^ij^ where Vij is the distance between particles i 
and j. Here Etrap = ^{(x^ - V^) cos{nt)Vrf/2r^ + {z^ - - y"^ /2)Uac/<P) is the electric 
field in a linear ion trap p , where z is along the trap axis. The radial and axial dimensions 
of the trap are tq and d. A necessary condition for stable trapping of noninteracting ions is 
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a g-parameter, qi = 2{Qi/mi)Vr:i/ {Q'^Tq) less than 0.9. However, it is well-known experimen- 
tally and theoretically that operation at significantly smaller q is favourable since rf-heating 



is less pronounced 
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25|. On the other hand, a lower limit is given by experimental con- 



siderations, gmin — 0.05. We therefore choose g-parameters in the range 0.05 - 0.4 for both 
atomic and molecular ions. As a consequence, the simultaneous storage of heavy molecules 
(msc ^ m\c) requires a molecular charge exceeding unity. 

Laser cooling of the Ic particles is described by the force Fl. In actual experiments, 
its strength is such that cooling may require minutes. Such durations correspond to 10^ 
or more rf periods, and cannot be simulated in high-precision MD when particle numbers 
are large. In order to compromise between a reasonable computing time and realistic laser 
cooling strength we have used stronger forces Fl to speed up cooling. Most simulations 
were performed with a simple linear viscous damping Fl = — /3x with friction coefficients in 
the range (3 = (1.2 — 8) ■ 10~^^kg/s. These are well below the maximum value (at optimum 
laser detuning from the cooling transition), /3max — tt^/I/A^ ~ 4- 10~^^kg/s for transitions at 
optical wavelengths A. For some simulations, we have used more realistic forces, see below. 

The diffusion of the Ic ion momenta in momentum space due to recoils upon spontaneous 
emission gives rise to the Doppler cooling limit. This diffusion is included in the simulations. 

The equations of motion are solved using a high-order Runge-Kutta method with adap- 
tive step size. Initial conditions were in part chosen so as to give initial temperatures below 
room temperature, again in order to reduce the computational effort. 

In order to characterize the state of the plasma, from the trajectories of the par- 
ticles we calculate, for each species, the average kinetic energy per particle, = 
(2A'^fj)~^ mj(vf (t)), where (■) is the time-average over one period of the rf field, the 
time-averaged secular energy per particle E^^^ = (2A''k)~^ (^^il^)) p5 the average 

interaction energy per particle (at the end of a rf period). Here k means Ic or sc, respectively, 
and YIk denotes summation over the corresponding species. The secular energy, where the 
micromotion oscillation is averaged out, can be taken as an indication of the temperature of 
the sample, since it arises from the "disordered" motion of the interacting ions in the time- 
averaged trapping potential. In contrast, the micromotion contribution to the total energy 
arises from a regular motion. We note that in simulations where the rf potential is replaced 
by the pseudopotential, i.e. the time-averaged trap potential V^pseudo(p) = p'^Q'^V^^/^AmQ'^rQ) 
experienced by the particles in radial direction, the cooling is moderately faster. This implies 
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that rf heating is nonneghgible even at the relatively low g-values chosen here. 
In the following we describe three mass regimes of sympathetic crystallization. 
Cooling of molecular hydrogen isotopomers One challenging goal in ultracold molecule 
studies will be the precision spectroscopy of the simplest (i.e. one-electron) molecules, the 
hydrogen ions H^, HD'^, D^. The experimental accuracy can potentially be improved by 
several orders, surpassing by far the current theoretical precision 23|. It will then become 
possible to test and challenge calculation methods, especially of the relativistic and QED 
contributions. Moreover, since the vibrational energies depend explicitly on the electron- 

their measurement might allow to determine the value of 



to-proton mass ratio nie/mp j21 
this fundamental constant by spectroscopic means, providing an alternative and potentially 
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more accurate approach than mass measurements in Penning traps 

A precise measurement will require cold trapped molecules in order to minimize Doppler 
broadening. The heteronuclear HD^ is of particular interest since it has dipole-allowed 
vibrational transitions 0] that could be excited by infrared lasers such as optical parametric 
oscillators or diode lasers. 

A simulation of sympathetic cooling of 5 HD~''-ions by 20 laser-cooled ^Be'''-ions is shown 
in Fig.l. The laser quickly cools the Be'''-ions to a liquid state, characterized by comparable 
Coulomb interaction energy and secular energy (plasma parameter F ~ 2). The atomic ion 
temperature remains constant while the molecular ions are sympathetically cooled. Only 
once the molecular secular energy becomes comparable to the atomic secular energy does 
the latter decrease further. The secular energy of the Be'''-ions finally reaches the Doppler 
limit. The secular energy of the molecules reaches that level significantly more slowly, since 
the cooling power of the atomic ions becomes smaller as they settle into the crystalline state. 



De understood by consider- 
32|. It is three times larger 



The spatial structure resulting from the cooling (Fig. ^p) can 
ing the different pseudopotentials felt by the two particle species 
for the lighter HD^-ions as compared to the ^Be'''-ions. The total energy of the ensemble 
is minimized if the HD'''-molecules lie on-axis. The Be'''-ions form a shell structure around 
them. It is moderately prolate, since the radial pseudopotential and axial potential are sim- 
ilar in strength (Be'^'-oscillation frequencies ujp/2'K = 340 kHz, Uz/^-k = 285kHz). The axial 
arrangement of the HD'''-ions is favourable for spectroscopic investigations, since on-axis 
the micromotion is zero, with a corresponding simplification of their transition spectrum. 
Since the Be'''-ions form a three-dimensional structure, their micromotion energy remains 
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relatively high compared to their secular energy (FigC^). This is because the off-axis lo- 
cations of the ions imply corresponding micromotion velocities, proportional to the radial 
distances. 

Simulations with a more realistic (i.e. weaker) cooling force js^ were also performed. 
Qualitatively, the same behaviour resulted, however, as expected, sympathetic crystallization 
was reached after a substantially longer time, about 10 times slower than in Fig^. 

A sample of molecular ions is not always sympathetically cooled in its entirety. For 
example, in a simulation of 40 Be'''-ions and 10 HD''"-ions, the final state of the system 
contained 6 on-axis crystallized HD'''-ions embedded in a prolate Be^-crystal, while 4 HD"*"- 
a) 
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FIG. 1: Simulation with Nic = 20 laser-cooled ^Be'''-ions and A'sc = 5 HD^ -molecules, a: Energies 
per particle, b: Spatial structure in the crystalline state (left: projection onto x-z-plane, right: 
x-y-projection); solid circles: ^Be+, circles: HD+. The trap parameters are qic = 0.13 Qsc = 0.39, 
rf frequency Q/ln = 8.5 MHz, 14f/^o ~ 17.6 V/mm^, U^c/d'^ = 30V/cm^. Viscous laser cooling 
with /3 = 2.4 X 10~22kg/s was used. 
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ions remained hot. It is clear that for small particle numbers the number of cooled particles 
depends on details of the initial conditions. 

The above simulations for HD"*" are of course also applicable to the astrophysically im- 
portant H;^. We have also performed simulations for other hydrogen ions. We have found 
(partial) sympathetic crystallization for all masses, from 2 amu (H^) to 5 amu (DT^). 

Cooling of dye molecules Dye molecules are interesting model systems for studies of com- 
plex (i.e. polyatomic) ultracold molecule manipulation because they are well-characterized 
and can easily be excited optically. Various dye molecules have masses exceeding a few times 
the mass of heavy atoms. We consider here the case of cooling Rhodamine 6G, a molecule of 
mass 493 that can be transferred to the gas phase singly charged by means of electrospray 
ionization. The coolant ion is chosen to be ^^^Ba"*". In this case, the larger mass-to-charge 
ratio of the molecules leads to a shell structure with the atomic ions in the center and on- 
axis. The small fraction of the molecules which is well embedded in the atomic ensemble 
cools and crystallizes on the same timescale as the atomic ions. The remainder experiences 
a much weaker cooling due to the absence of a "cageing" effect. The timescale for cooling 
and crystallization of all molecules was found to be 10 times larger than for the atomic ions. 

Cooling of large molecules Large molecules such as amino acids or proteins can be trans- 
ferred into vacuum in high charge states by electrospray ionization 
protons are attached to molecules that emerge from a solution spray. Large molecules are 
produced with a distribution of charge states, with typical specific charges in the range 
m/Q = 700 — 1500. For example, the charge states for a protein of mass 17,000 can range 
from 11 to 25 2^. Such specific charges are sufficiently near those of laser-coolable heavy 
atomic ions (^^'^Ba'*', ^^^Yb"^, ^^^Hg"^). Since the g-parameters ~ Qi/rrii are within an or- 
der of magnitude, simultaneous trapping is possible. Moreover, the large charge ensures a 
good coupling between molecular ensemble and atomic ensemble. Finally, the pseudopoten- 
tial of the molecular ions (~ Qg^/rrisc) can be comparable or even larger than that of the 
atoms, thanks to the large molecular charge. A good spatial overlap of the subensembles 
is therefore ensured. FigEl shows two simulations for molecular ions of mass 20,000 amu 
and charge 20 e. As can be seen, the molecular ions are strongly cooled as far as their sec- 
ular energy is concerned. However, a large micromotion energy remains in the crystallized 
state because of the nonaxial arrangement of the molecular ions and their large mass. This 
nonaxial arrangement occurs in spite of the fact that the pseudopotential of the molecules 



2J|. In this method. 
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FIG. 2: Sympathetic cooling of 10 molecular ions of mass 20000 and charge 20 e by 30 '^^'^Ba+- 
ions. (a) Energies, (b) Spatial structure, (c) Spatial structure for a simulation where the axial 
potential is weaker. Parameters: qic = 0.33 Qsc = 0.045, r2/27r = 1.6 MHz, Vr{/2rQ = 23.8V/mm^, 
Udc/d? = 12V/cm^ for (b), 3V/cm^ for (c). Viscous laser cooling with /? = 4.8 x 10~^^kg/s was 
used. 

is larger than that of the atoms. Indeed, the Coulomb interaction energy is an important 
contribution that influences the overall arrangement of the particles. This can be seen by 
noting that the Coulomb energy of a given spatial structure is not invariant under particle 
exchange Ic^sc, due to the differing molecular and atomic charges. For a sufficiently large 
axial potential (FigI2t,b) the total energy is minimized by forcing the molecular ions off axis 
resulting in a shell structure. For smaller axial potential (FigI2fc) the structure is string-like, 
but molecular ions that are adjacent bulge out of the axis because of their strong repulsion. 
Where several molecular ions are crystallized in adjacent positions, they form sections of 
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a helix 3^. Here again we find that molecules that are well embedded in the atomic ion 
ensemble, i.e. those that are individually located between atomic ions (e.g. the two sc ions 
in Fig|2t furthest from the center in z-direction) , cool as fast as the Ic ions while the off-axis 
molecular ions cool much more slowly. 

We have also performed simulations for molecular masses 2000, 5000, and 10,000 amu, 
with the same mass to charge ratio of 1000 amu/e. We find a similar behaviour of crystal- 
lization as in the examples above. As the charge becomes smaller only molecular Coulomb 
structures of the shell-type occur, and due to the weaker Ic-sc interaction, a fraction of the 
molecules remain uncooled. 

The above examples show that sympathetic crystallization is possible for ion masses 
significantly larger or smaller than the coolant ion mass. Based on the present and previous 
j23[ results, it can be stated that simply charged molecular ions of mass between that of 
^Be^ and ^'^^Ba'*' can also be sympathetically crystallized. 

In conclusion, we have shown that a wide variety of charged molecules can be trapped and 
cooled to niK temperatures with masses from 2 to 20,000 amu, i.e. from diatomic molecules 
to polymers and proteins. The molecules are incorporated into an ordered structure. The 
spatial arrangement within the Coulomb crystal depends on the masses and charges of the 
coolant and the molecular ions. We have pointed out that sympathetic crystallization should 
be very advantageous for precision spectroscopy of the various ions of molecular hydrogen. 

Since emphasis was placed on using a laser cooling force of realistic magnitude and on 
employing a high-precision numerical code, the simulations were performed for small particle 
numbers (up to 50). Experimentally, this regime is both accessible and suited for detailed 
studies of molecular properties. In future, the simulations can be extended to larger numbers, 
by developing faster but approximate algorithms or using more powerful computers, and to 



othe. types of t..aps (.u..p„.e.rf t.aps B and Pe,.,„g t.aps), TMs ..„„M allow detaUed 

comparisons of dynamics and structure of sympathetic crystallization between theory and 
experiments currently under way. 
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